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SUMMARY 


A simple unified numerical method applicable to non-planar subsonic air- 
craft has been developed for calculating either the induced drag for an arbi- 
trary loading or the optimum aircraft loading which results in minimum induced 
drag. The method utilizes a vortex lattice representation of the aircraft 
lifting surfaces coupled with the classic equations and theorems for computing 
and minimizing induced drag. Correlation of results from the numerical method 
with non-planar solutions obtained from other more complex theories indicates 
very good agreement. Comparison of the induced-drag computations using the 
numerical method with experimental data for planar and non-planar configura- 
tions was also very good. 


INTRODUCTION 


Over the past few years, increased attention has been focused on improving 
aircraft performance. One method to improve performance is to lower the air- 
craft induced drag. This can be accomplished by more efficient design of con- 
ventional configurations or by developing new and unique designs whose intent 
is to minimize induced drag. Typical of new configurations that have been 
developed for this purpose are the Lockheed boxplane and the Whitcomb winglet 
configuration. 

The aircraft lifting surfaces for conventional aircraft as well as for new 
configurations are generally non-planar in design. To achieve a minimum in- 
duced djag, these non-planar surfaces must be designed to support the required 
optimum loads as specified by classical theory (refs. 1 and 2). Unfortunately, 
the use of classical theory to determine the design loads is quite cumbersome 
since rather complex conformal transformations must be utilized. Thus, a 
simple inexpensive method is required to determine vrhat the ’’design to” loading 
of a no.*“ planar configuration should be to minimize the aircraft induced drag. 
Furthermore, for conditions where the aircraft is not operating at design con- 
ditions, an analysis method is required to quickly assess the magnitude of the 
aircraft off-design induced drag. Also, methods of this type are of particular 
importance in making configurational trade-offs. 

The objective of this paper is to present a simple unified numerical 
method applicable to subsonic non-planar aircraft for the rapid calculation of: 
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1. the Induced drag for an arbitrary aircraft loading or 

2. the optimum aircraft loading which results In minimum Induced drag. 

The paper will Include a discussion of the fundamental theoretical concepts on 
which the method Is based, followed by the theoretical formulation of the nu- 
merical calculation procedure. Computations will be made using the method and 
will be compared to existing theoretical solutions and to experimental data. 
This will be followed by an Illustration of the utility of the method for 
making configurational trade-offs by comparing the loading and Induced drag 
results for various types of wing additions such as wlnglets or wlng-tlp 
extensions. 
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geometric Influence function 
aspect ratio, b^/S 
reference span 

local chord of lifting surface 
average chord (S/b) 

section load coefficient normal to load perimeter 

lift coefficient 

Induced drag coefficient 

bending-moment coefficient 

Induced drag (D^ = Cj)^ qS) 

efficiency factor 

resultant force of lifting surface 

length of wing addition 

unit normal vector parallel to Z axis 

length of load perimeter 

lift force (L = Cl qS) 

Mach number 


50 




I 




number of lifting elements on load perimeter 
unit vector normal to load perimeter 
section load normal to load perimeter (N * q Cjj c) 
free-stream dynamic pressure 
semi-width of vortex pair 

nondimens ional semi-width of vortex pair (s»2s'/b) 


reference area 


side force 


induced velocities 


resultant induced velocity 
velocity normal to load perimeter 
free-stream velocity 
lifting element coordinate system 
aircraft coordinate system 
aircraft spanwlse center of pressure 
circulation (eq. (10)) 
rotation angle in the Y-Z plane 


Superscript: 


Subscripts: 


indicates vector quantity 


number designating a vortex pair that model a particular lifting 
element 

number designating a control point on a particular lifting element 


winglet 
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BASIC THEORETICAL CONCEPTS 


Fundamental to the development of the present model is the representation 
of the aircraft non->planar lifting surfaces by a system of rectangular horse- 
shoe vortices (ref. 3). The Induced drag for a given loading or the optimum 
loading for minimum Induced drag can be calculated for any arbitrary non-planar 
aircraft at subsonic speeds utilizing this vortex representation and the 
following basic law and theorems: Hunk's Theorems I to III« Blot-Savart Law, 

and the Kutta-Joukowskl Theorem. 

Hunk's fltst theorem (ref. 1) can be stated as follows: 

The total induced drag of any multiplane system is unaltered 
if any of the lifting elements are moved in the direction of 
motion provided that the attitude of the elements is adjusted 
to maintain the same distribution of lif among them. 

This theorem Is commonly referred to as Hunk's stagger theorem. An Illustra- 
tion of this theorem Is shown In figure 1. Several practical applications can 
be deduced from this theorem. First, the chordwise distribution of pressure 
does not affect the theoretical Induced drag of the aircraft If constant 
section lift Is maintained. Second, wing sweep does not effect the theoretical 
Induced drag as long as the spanwlse distribution of lift Is constant. A third 
application Is that the load from a system of multi-surfaces (l.e. wing and 
horizontal tall) with the same projection in the Y-Z plane can be made equiva- 
lent to a single surface for the purpose of calculating Induced drag. 

In the following theoretical development, use will be made of Hunk's first 
theorem to lump the chordwise distribution of vorticity into a single chordwise 
load and to translate all loads into the 0,Y,Z plane (fig. 1). 

Hunk's second theorem (ref. 1) is illustrated in figure 2 and can be 
stated as: 


In calculating the total induced drag of a lifting system^ 
once all the forces have been concentrated into the plane 
OiY,Zf we may, instead of using the actual values of the 
velocity normal to the lifting elements [VYi(x,y,z)] at the 
original points of application of the forces, use one-half 
the limiting value of the normal velocity [Vfi(‘°,y,z)] for 
the corresponding values at points P(0,y,z). 

This theorem allows the computations to be done in the Trefftz plane (down- 
stream Infinity) rather than in the real plane. In the subsequent theoretical 
derivation, this fact will be utilized to make all the computations in the 
Trefftz plane, thereby greatly simplifying the calculations. 

The third theorem given by Hunk (ref. 1) is presented as follows: 


52 



When all the elements of a lifting system have been trans- 
lated longitudinally to a single plane ^ the induced drag 
will be a minimum when the component of the induced 
velocity normal to the lifting element at each point is 
proportional to the cosine of the angle of inclination of 
the lifting element at that point. 

This theorem is illustrated in figure 3 and can be summarized in equation form 
as: 


Vn * Wq cose (1) 

For a horizontal lifting element it can be seen from equation (1) that the 
normal velocity (downwash) across the span is equal to a constant (fig, 3). For 
a vertical plane (0=90°), the normal velocity (sidewash) must be equal to zero 
for minimum induced drag. The physical interpretation of this theorem will be 
further illustrated in a subsequent section. 

Equation (1) will be utilized in the following theoretical development as 
the boundary condition necessary to achieve a minimum induced drag and hence an 
optimum aircraft loading. 

The basic equation for calculating the aircraft-induced drag can be 
derived by applying the Kutta-Joukowski theorem in the drag direction. By 
virtue of Munk’s theorems, the calculations can be accomplished in the Trefftz 
plane rather than the real plane. Thus, the equation for induced drag ex- 
pressed in terms of the Trefftz plane variables and using vector notation is: 


2V 


(p V • Nn dil 


( 2 ) 


Equation (2) along with the induced velocities in the Trefftz plane de- 
rived from the vortex model of the lifting surfaces w:Ill comprise the basis for 
the induced drag computation. 


PHYSICAL INTERPRETATION OF THEORETICAL CONCEPTS 


To provide a better physical understanding of the computation of induced 
drag and the calculation of the optimum loading for minimum induced drag, the 
theoretical concepts discussed in the previous section will be illustrated 
using a wing-winglet configuration. In figure 4, the sources of induced drag 
for a wing-winglet combination are shown. These are: 

o Drag due to the induced flow by the wings on the wing 

o Drag due to the induced flow by the wings on the winglet 
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o Drag due to the Induced flow by the wlnglets on the wlnglet 

o Drag due to the Induced flow by the wlnglets on the wing 

For simplicity, the effects of symmetry are Included in the sources of Induced 
drag shown and are not delineated separately. 

In figure 4(a), the effect of the wing Induced flow is shown. The wing 
under positlve_load produces a downwash on Itself which results in the wing 
force vector, F, tilting rearward by an angle oi. The wing force vector, P, is 
perpendicular to the resultant, V, (Kutta-Joukowski theorem). The rearward 
rotation of the force vector results in a wing-induced drag. A sldewash is 
also produced by t «e wing at the wlnglet location. As can be seen in figure 
4(b), the sldewash from the wing combined with the free-stream velocity pro- 
duces a tilt forward of the wlnglet force resulting in a thrust component. 

In figure 4(c), the induced drag resulting from the sldewash of the win'^- 
let on itself is presented. This results in a rearward tilting of the wing'.et 
force vector and an attendant induced drag. It should be noted that the d:*.rec- 
tlon of the wlnglet force vector is consistent with a positive (upload) or the 
wing. The wlnglet also Induces an upwash on the wing. In figure 4(d), it can 
be seen that this upwash rotates the wing force vector forward producing a 
thrust force. 

The results from figure 4 are summarized in figure 5, where all the in- 
duced velocities are combined. For minimum Induced drag, equation (1) indi- 
cates that the velocity normal to the wlnglet must be equal to zero (6 ° 90°) . 
This can be seen to occur when the sldewash produced on the wlnglet by the wing 
exactly cancels the sldewash produced by the wlnglet on itself. In other 
words, the induced angle of attack (a^) of the wlnglet is zero. The induced 
drag of the wing is also minimized by the presence of a wlnglet since the wlng- 
let causes a reduction in the net downwash at the wing; and, hence, the induced 
angle of attack is reduced. 


DERIVATION OF NUMERICAL METHODS 


Vortex Model 

By virtue of Hunk's theorems, the calculations for induced drag and the 
optimum loading can be accomplished in the Trefftz plane. This fact consider- 
ably simplifies the calculation problem since the method will not be a function 
of the longitudinal coordinate. The projection of the aircraft non-planar 
lifting surfaces i the Trefftz plane will be referred to as the load 
perimeter. 

In the real plane, the aircraft lifting surfaces will be represented by a 
system of horseshoe vortices. The equations describing the Induced velocities 
in the Trefftz plane at a control point P(«>,yi^,Z£) (fig. 6) due to a horseshoe 
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vortex located In the real plane at a point P(Xj,yj,Zj) are given below as de- 
rived from the Biot-Savart Law (ref. 3): 


v„ 2 uv^\ri R2I 

, 1 f (y' -s') (y‘ +s') \ 

27t \ Ri R2 / 

where 

Rl = (z')2 + (y' -s')2 (6) 

R 2 = (z')2 + (y' +s')2 (7) 

y' = (Yi-yj) cos0j + (z^-Zj) sin6j (8) 

z' = “ (yi~yj) sinOj + (z£ - zj) cos0j (9) 

and the circulation by virtue of the Kutta-Joukowski theorem is given as 

(cn c)-j 

( 10 ) 

00 

Inspection of equations (3) to (5) indicates that there is no contribution 
from the horseshoe bound leg in the Trefftz plane and the induced velocities 
are not dependent on x. The resulting model then reduces to describing the 
load perimeter in the Trefftz plane (fig. 6) by lifting elements that are 
represented by a trailing vortex pair having a circulation of equal magnitude 
but of opposite rotation. For each lifting element there is an associated 
control point located midway between the pair of vortices. 

Induced Drag Calculation 

The basic equation for calculating the induced drag for an arbitrary non- 
planar lifting system was given in equation (2) as: 

D^ = -::h- I V • Nn d£ 


The integral is a circuit integral taken around the perimeter of the pro- 
jection of the lifting system in the Trefftz plane. The vector V is the 
resultant induced velocity vector in the Trefftz plane from all vortices on the 
load perimeter. The vector n is a unit vector, normal to the load perimeter. 






Reducing equation (2) to coefficient form, nondlmensionalizlng the lifting 
element length (Z) by the reference semispan, and using the relationship for 
the average chord, the following result is obtained 


Cdi 




( 11 ) 


Writing the above in the form of a svm and assuming S5rmmetry about the X-Z 
plane 




1 y ^^n c)j 

2 iii cav 


[V cos(v,a)]. 

^ A(2)l/b)i 

00 


( 12 ) 


where m equals the number of elements that comprise the load perimeter. 

Writing equation (12) in terms of the nondimens ional lifting element semi- 
width (s) and noting thac 

Vj^ = Vcos(V,n) (13) 


the expression for induced drag can be written as 


^Di 




(14) 


The velocities normal to the lifting elements (Vj^) can be determined by 
utilizing the expressions for the induced velocities in equations (A) and (5). 

From the geometry of figure 6, the normal velocity at due to a vortex 
pair at Pj can be expressed in terms of the induced velocities as: 

w± Vi 

— = — cos (0^ - 0j) - — sin(0i - 0j) (15) 


Combining equations (4), (5), (10), and (15) yields the expression for the 
total normal velocity at the control point Pi due to vortices at all points P j : 


^co j=i <^AV 


^AV / (v‘ - s') 
4ti \ Ri 


(y' +s') \ 

R2 / 


cos(6j^ ” 


+ 


^AV 

4ti \Ri 



sin(0£ - 0j) 


(16) 


The portion contained in curvy brackets is only a function of the projected 
aircraft geometry in the Trefftz plane and will be denoted by . Thus, in 
terms of the geometric influence function A^ j , 
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Vhj ro (Cn c) , 

^ = I ^ A.. (17) 

V«> c^V ^ 

Substitution of equation (17) Into equation (14) yields the final expres- 
sion for the induced drag: 


r r 


mm/ 

)i= ^ I 

1=1 i=l ' 


•) (siHAij) 


The independent parameters in equation (18) are the loadings normal to the 
load perimeter, the lifting element semi-widths, and the geometric Influence 
function. The loading normal to the load perimeter will be considered input to 
the present method. The norma] loading can be determined from any available 
non-planar lifting surface calculation procedure such as in reference 3 or from 
experimental data. The lifting element semi-width is also considered as input. 
The geometric influence function as has been previously mentioned, is a 

function of the input aircraft geometry. 

For an arbitrary applied load, the lift can be determined from the follow- 
ing expression 


L = 0 Nn . k dii 


Expressing equation (19) in coefficient form and writing as a sum 


Cl = 2 


i (^) 

j = l ' ‘^AV / 


(Sj) cos(iij, k) 


Since 


cos (nj , k) = cos 0 j 

the final expression for the lift coefficient is given by 
m / (cn c) \ 

^ ®J 


The bending-moment coefficient at the X-axis can be expressed as 


2 >1 ^ <=AV V'' 


Yi z- 

b72 j b7 


^ sin0jj 


The spanwise center-of-pressure location can be determined from the following 
equation : 


I 
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ycp 
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The aircraft efficiency factor can be calculated from the following 
standard equation: 


IT AR Cr 


Optimum Load Calculation 

The expression for the total velocity normal to a lifting element was 
given in equation (17) as: 

Vni m (cp 

j-i *^AV 

According to Hunk's theorem III, the loading for minimum Induced drag Is 
obtained when the distribution of normal velocity satisfies equation (1): 


V T 


where w^ Is a constant. Using equation (1) as a boundary condition and combin- 
ing It with equation (17), there results: 


Wo “ (<^n c) 

7" - I 

» j-1 *^AV 


where the loading In equation (2b) Is the optimum loading. Using square 
brackets to Indicate matrix notation, equation (26) can be written as: 


vJV 
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!!^ 


Solving for the optimum loading 


(cn c) 


^ A 
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The value of the arbitrary constant, w^, can be determined from equation 
(22) by specifying the aircraft lift coefficient. 

In summary, to determine the loading for minimum induced drag, only the 
lift coefficient and aircraft geometry are required for Input. Once the 
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locidlng has been determined, the minimum Induced drag, the spanwise center-of- 
pressure location, and efficiency can be determined in the manner previously 
presented (eqs. (18), ( 24 j, and (25)). 


CORRELATION OF METHOD 


The theory described in the previous section for calculating the aircraft 
Induced drag for a given loading or th.? loading for a minimum Induced drag has 
been coded for use on Lockheed computers. In this section, computations using 
the present method will be compared to other theoretical solutions and to 
experimental results. 


Induced- Drag Correlations 

Numerical solutions for the aircraft efficiency for a monoplane of aspect 
ratio eight are shown In figure 7 for various values of lifting element widths. 
For this example, the widths of the elements over the load perimeter were held 
constant. A more efficient result could have been obtained If, for Instance, a 
cosine spacing of the elements had been used. The Input loading on the mono- 
plane was specified to be elliptical. The exact solution for the efficiency 
factor (ref. 1) on an elllptlcally loadvsd monoplane Is, of course, 1.0. As can 
be seen, the numerical solution approaches the exact value as the width of the 
elements become smaller. For a lifting element width equal to .01 (b/2) , the 
error was approximately 0.5% In efficiency. 

A similar calculation was made for an aspect ratio eight biplane with 
wings of equal span and a helght-to-span ratio of 0.5. The biplane was loaded 
optimally utilizing the loadings derived In reference 1, based on transforma- 
tion theory. The numerical calculation for efficiency factor was 1.6307 (using 
constant elements of 0.0125(b/2) in width) compared to the value of 1.6260 
given In reference 1. The resulting difference was approximately 0.3%. 

In figure 8, Induced-drag results calculated using the present method are 
compared to the experimental results for an advanced Lockheed transport de- 
signed to cruise at 0.95 Mach number and at a lift coefficient of 0.47. The 
spanwise loading for the aircraft was obtained from an available lifting 
surface program similar to that In reference 3. As can be seen, the agreement 
between theory and experiment Is very good over a range near the design lift 
coefficient. 

A further correlation example is presented in figure 9 where numerical re- 
sults are compared with experimental data for a non-planar Lockheed boxplane 
configuration. Again, the loading was obtained from lifting-surface theory. 

As indicated, good agreement Is obtained. 


Optimum Load Correlations 

In figure 10 the optimum loading Is presented for an aspect ratio eight 
monoplane as calculated from equation (28) using constant elements of .01 (b/2) 
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In width. Also shown Is the classic optimum result for a monoplane an 
elliptical loading. The correlation can be seen to be very good. The error in 
Induced drag of the computed result was approximately 0.5%» which is consistent 
with the results in figure 7. 

In reference 2, an optimum loading for a wing with a winglet is presenLed 
as derived from transformation theory. This solution is compared to the result 
calculated using the present method in figure 11. Good agreement between the 
two methods is obtained. 


ILLUSTRATIVE USE OF METHOD 


Taken together, the present numerical method provides a unique tool for 
understanding the sources of Induced drag and making configuration trade-offs 
to achieve an overall aerodynamic as well as structurally optimum aircraft. An 
illustration of using the method to provide additional understanding into the 
basic sources and mechanisms of Induced drag is presented in figure 12, where 
the magnitude of the induced-drag components for a wing-wlnglet configuration 
(fig. 4) with optimum load is presented. From the figure, it can be seen that 
the induced-drag contribution from the wing on the winglet and :he winglet on 
the winglet are of the same magnitude and cancel each otner. This is, of 
course, the result previously illustrated in figure 5. 

As a result of design or structural constraints, the aircraft may not be 
able to achieve the op 'muir loading for minimum induced drag. The penalties 
that incur from the use of non-optimum loadings can be quickly assessed using 
the present method. This is illustrated in f.’.gure 13, where the winglet load- 
ing for a wing-winglet configuration is varied. As can be seen the Induced 
drag is sensitive to certain types of changes (non-optimum 1) where it is not 
to others (non-optimum 2). 

The present method can also be used to quickly make configuration trade- 
off3. This is illustrated in figures 14 and 15. In figure 14 the paiametric 
effect of wing additions on induced drag and on the wing spanwise center-of- 
pressure location as calculated from the present metnod are presented. If, for 
instance, it was desired to find a configuration which would give the maximum 
induced drag reduction for a minimum outboard shift in wing center of pressure, 
this can quickly be determined by replotting the parametric data of figure 14 
in the form of figure 15 and the result determined. 


CONCLUDING REMARKS 


A unified numerical method applicable to non-planar subsonic aircraft hcs 
been developed for the purpose of calculating the induced drag for an arbitrary 
loading or the optimum aircraft loading that gives minimum induced drag. 
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Use of the numerical method has indi::ated that: 


(1) the method is simple and easy to use 

(2) induced drag and optimum loading results from the numerical 
method correlate very well with non-plauar solutions obtained 
from more complex theories 

(3) numerical induced-drag predictions ate in good agreement with 
experimental data for planar and non-planar configurations 

(4) the numerical method provides both analysis and design capa- 
bility which allows the designer to make rapid configuration 
assessments and trade-offs for the purpose of achieving an 
overall aerodynamic as well as structurally optimum aircraft. 
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Figure 4.- Sources of induced drag for a wing/winglet configuration 
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Figure 6.- Induced velocities in the 
Trefftz plane. 
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AR =8.0 



Figure 7.- Variation of efficiency for optimally loaded monoplane 
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Figure 9.- Induced drag for a Lockheed boxplane- 
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Figure 10.- Comparison of numerical and exact optimum 
span loadings for a monoplane. 




Figure 11.- Comparison of theoretical optimum loadings 
for a wlng/winglet configuration. 
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Figure 12.- Illustration of induced drag calculations for 
a wing/winglet configuration. 
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Figure 13.- Effect of changes in winglet loading on the ind ^ed drag 
of a wing/winglet configuration. 




WING CENTER-OF-PRESSURE 
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Figure 14.- Basic theoretical effects due to wing additions. 
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LIFT = CONSTANT 
MINIMUM DRAG LOADING 



WING CENTER-OF-PRESSURE CHANGE 
(°/c SEMI SPAN) 


Figure 15.- Theoretical effect of wing additions 
for various dihedral angles. 




